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Abstract 



Possible Lorentz invariance violation (LIV) has been investigated for a long time based on observations of GRBs. These 
^1 .arguments relied on the assumption that photons with different energy are emitted at the same place and time. In 
|>^this work, we try to take account of the intrinsic time delay A£j nt between emissions of low and high energy photons 
,by using the magnetic jet model. The possible LIV effects are discussed in a unified scenario both for long and short 
^> "Fermi-detected GRBs. This leads to a unique quantum gravity energy scale M\c 2 ~ 1.0 x 10 20 GeV respecting the linear 
dispersion relation. 
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1. Introduction 



Lorentz invariance is one of the most important corner- 
q -stones of modern physics. Recently, the OPERA collab- 
*H [oration reported that the GeV neutrinos propagate faster 
t/3 than the speed of light [l[. Although other independent 
i ^ / experiments should be performed to verify the superlu- 
minal phenomenon, it is valuable to ask whether Lorentz 
invariance violation (LIV) happens in high energy scale. 



A favored way to test LIV is to study the most explosive 



events in the present universe: gamma-ray bursts (GRBs). 
t^J- -The Fermi satellite has observed several GRBs with phlo- 
em "ton energy > 100 MeV in recent years. Fermi carries two 
* .instruments: the Gamma-Ray Burst Monitor (GBM) and 
<^ [the Large Area Telescope (LAT) , which detect the energy 
(N -band 8KeV - 40MeV and 30MeV - 300GeV, respectively. 
T j~^ [An interesting feature of the observation is that GeVpho- 
;> tons arrive several seconds later than MeV photons [2t5|- 
| One possible explanation is given by quantum gravity 
.effects. Some quantum gravity theories predict that high 
C$ energy photons may interact with the foamy structure of 
the space-time, thus photons with different energy prop- 
agate with different velocities 6-9]. Such effects can be 
accumulated after photons travel a cosmological distance. 
In these theories, high energy photons are subluminal. But 
there are still some other theories which show that high en- 
ergy photons can be superluminal. For example, LIV can 
also be induced from the geometry of the space-time itself, 
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such as the Finsler geometry 1(J- 13] ■ One can expect that 
the velocity of photon also depends on the energy, and may 
be superluminal. 

A straight forward way to test LIV is studying the in- 
dividual GRB, several papers have discussed the upper 
limits of variations of the light speed 13j, [14J. Ellis et 
al. have proposed a data fitting procedure to test LIV 
effects [tj, Il5l - [r7| . The linear fitting function is expressed 
as At bs/(1 + z) = a\ J \YK{z) + b, where K(z) is a non- 
linear function of the redshift measuring the cosmological 
distance, and b represents the ignorance of the intrinsic 
time lag. A statistical ensemble of GRBs was used to fit 
values of oliv and b, and no strong evidence of LIV was 
found. The explicit form of K(z) depends on the cos- 
mological model. Biesiada and Piorkowska applied this 
procedure to various cosmological models [18(. Shao et 
al. used this method to discuss four Fermi-detected GRBs 
171 ] , which we will discuss in the present work. Neverthe- 



less, all these investigations concentrated on the time lag 
induced by LIV, the intrinsic time lag which depends on 
the emission mechanism of GRBs was not considered. 

On the other hand, without considering LIV effects, sev- 
eral mainstream GRB models were proposed to explain 
the delayed arrival of GeV photons ljj-22|. Meszaros 
and Rees presented a magnetic dominated jet model to 
explain this phenomenon. The MeV photons can escape 
the plasma when their optical depth decreases to unity at 
the photosphere radius. While the GeV photons are pro- 
duced by the nuclear collision between protons and neu- 
trons, which happens at a large radius (compared to the 
photosphere radius) [20| . Duran and Kumar considered 
that photons are emitted by electrons via the synchrotron 
radiation, it consumes more time for electrons to be accel- 
erated in order to radiate GeV photons [2l[. Bosnjak and 
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Kumar proposed the magnetic jet model, the time delay 
depends linearly on the distance where the jet is launched 
22]. 

In this work, we argue that the observed time lag for 
two photons with energy -Ehigh and E\ ow consists of two 
parts, 



Aiobs = A^liv + Af i. 



(1) 



where At; nt denotes the intrinsic emission time delay, and 
AtLiv represents the flying time difference caused by LIV 
effects. In section^ we take use of the magnetic jet model 
in Ref. 22j to estimate Ai; nt . The LIV induced time lag 
A^liv is given, and the quantum gravity energy scale is 
discussed in Sec. [3l Finally, the discussion and conclusions 
are given in Sec. 2] 

2. Magnetic jet model 

In the magnetic jet model, photons with energy less than 
10 MeV can escape when the jet radius is beyond the 
Thomson photosphere radius, i.e., the optical depth for 
low energy photons is tx ~ 1. However, GeV photons will 
be converted to electron-positron pairs at this radius, and 
can escape later when the pair-production optical depth 
r 77 (£') drops below unity. 

The bulk Lorentz factor of an expanding spherical fire- 
ball increases with the radius roughly as F oc r, until 
reaching a saturate radius r s where the Lorentz factor is 
saturated |23| . However, for an expanding jet with small 
ejecting angle, the effective dynamical dimension is one. 
The bulk Lorentz factor increases with the radius roughly 
as 



22] 




for ro < r < r s 
for r > r., 



(2) 



where ro w 10 7 cm is the base of the outflow, which rep- 
resents the distance from the central engine where the jet 
is launched. 77 is the final bulk Lorentz factor of the jet. 

The optical depth for the photon-electron scattering is 
defined as 

f°° dr' 

tt(t) = J ^a T nT, (3) 

where ctt = e 4 /(67T£gC 4 mg) » 6.65 x 10~ 25 cm 2 is the 
Thomson scattering cross-section. The baryon number 
density in the observer frame is n ~ L/47rr 2 TO p Fc 3 <7o, 
where L is the isotropic luminosity, m p is the mass of 
protons, (To = L(ro)[l + a(ro)], and <j{ro) is the initial 
ratio of the magnetic and baryon energy densities. When 
t t{t) = 1) the Thomson photon sphere radius is 



LtL 

ro 



1.36 x 10 5 L^ 5 <7o 3 /5 r, 



3/5 
0,7 ' 



(4) 



where L 52 = i/10 5 erg • s , oo,3 = fo/10 , and r 0i 7 = 
ro/10 7 cm. We use cgs units for numerical values here and 
after. 
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Table 1: The observed parameters of four Fermi-detected 
GRBs. E p is the energy at the peak of vf v spectrum. £ , i os ,54 is 
the isotropic equivalent energy in unit of 10 54 ergs. T90 is the 
GRB duration which 90% of the counts are above background. 
2 is the GRB redshift El. 



The optical depth for a photon of energy Eq to be con- 
verted to while traveling through the jet at a radius r, 
is given by [22} 



r±(E ,r) 







(3-1 ) 471TT 2 (1 + z^-^EpC 



E P E 



0-1 



(5) 

where <j 11 = 6 x 10 26 cm 2 is the cross-section for photons 
producing just above the threshold energy, E p is the 
peak energy of the vF v spectrum, L >p is the frequency- 
integrated luminosity above E p , and j3 ~ 2.2 is the photon 
index of the spectrum above E p . Setting t± — 1, the pair- 
production photonsphcre radius is given by 
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(6) 

where E p - 6 = E p /MeV, and £ -4 = So/100MeV. 

In the observer frame, the relative time delay between 
MeV photons and GeV photons equates to the time for 
the jet to propagate from r p to r 77 [22j 



At 



3r Q (l + z) 
2c 



<{Eo) 



'•0 



1/3 



l/3-i 



(7) 



Taking use of Eq.© together with Eq.© and Eq.©, one 
can calculate the time delay for the arrival of GeV and 100 
MeV photons relative to MeV photons. Thus, the intrinsic 
time delay is At; nt = Ai(-Ehigli) — Ai(£i ow ). The observed 
values of E p , £ , i os ,54, T90 and z for four i^ermi-detected 
GRBs are taken from [22J , and are listed in Table [TJ 



3. Test of LIV effects 

As is mentioned above, some quantum gravity theories 
predict that photons with different wavelength propagate 
0, ES HI- The non-trivial space- 



in different speed 
time structure may affect the propagation of photons, so 
high energy photons may arrive later than low energy ones. 
Many works have studied LIV effects of the high energy 
photons in this direction 0, [3, 2^- 28| . 



2 



Consider two photons emitted at the same time and 
place, within the LIV phenomenology, the arrival time de- 
lay between them can be written as 25J, [27 1 



_l + n ( AE n 



(8) 



where n = 1 or n = 2 denotes the linear or quadratic 
correction to the dispersion relation, and D n is defined to 
be HSUlI 



(1 + z') n dz' 



Ho Jo y/riMii + z') 3 + n A 7 



(9) 



where Ho — 72 km sec -1 Mpc" 1 is the Hubble constant, 
fijvf and Q\ are the present values of the matter den- 
sity and cosmological constant density, respectively. In 
the standard cosmological model, (^m^a) are given by 
observations as (0.3,0.7) 27]. For n = 1, the time de- 
lay depends linearly on the variation of energy, which we 
will consider in the follow. In this case, the effective LIV 
energy scale is 



M lC 2 = 



, AED 1 



cAt 



LIV 



(10) 



In Ref.[l7J, Shao et al. took use of 4 GRBs the same as 
in Table [1] to predict LIV effects, the observed time delay 
Ai bs/ (1 + z) v.s. K(z) was plotted, where K(z) is defined 
as 



K(z) = 



AE D x 
(l + z)~' 



(11) 



In their plot, the three long bursts were found to be near 
one line. However, the short burst was not fitted well by 
the same line. The intercept of the line was interpreted as 
Ati n t/(l + z), and the slope of the line can be interpreted 
as 1/Mic 2 . Ati nt /(1 + z) was found to be negative, which 
means that high energy photons are emitted earlier than 
low energy ones. This conflicts with standard GRB mod- 
els. Taking into account of the four GRBs, the quantum 
gravity energy scale is estimated to be M\ c 2 ~ 2 x 10 17 
GeV for the linear energy dependence. In their work, the 
intrinsic time delay is assumed to be the same for three 
long bursts. Besides, their result strongly depends on the 
artificial choices. For example, if the 33.4 GeV photon is 
replaced by the 11.16 GeV photon in GRB 090902b, the 
three long bursts can not be fitted well by one line. 

In the follow, we first calculate Aii nt by using Eq.©. 
Then, combining with the observation data, we give A^liv 
by the relation given in Eq.(H]). In Ref.[25|, the authors 
used a statistical method to determine the observed time 
difference Ai b s between photons with energy -Eiow and 
-Eiiigh- In our present work, i?i ow is taken to be 100 MeV, 
since the onset of 100 MeV photons can be read directly 
from the data of the LAT monitor. -Ehigh is the energy of 
the most energetic photon in each GRB. One exception is 
that the second energetic photon with E^igh = 11.16 GeV 
in GRB 090902b is chosen. The most energetic 33.4 GeV 



photon arriving at 82 s is excluded due to its incoincidence 
with the main burst. In Ref. 17J, Shao et al. selected this 
event to estimate LIV effects without considering the cen- 
tral engines and emission mechanism. In the magnetic jet 
model, if we believe that the delayed GeV photons are due 
to the high optical depth, the arrival time of the 33.4 GeV 
photon should be 3 s later than the 11.16 GeV photon. 
However, the observed time interval 70 s is far beyond the 
model's prediction. This photon may be due to the in- 
elastic collision of protons and neutrons [l9L I20I ] , and it is 
quite possible that this individual event happens when the 
jet encounter the interstellar medium. 

With Eq.0, we can estimate At int ~ 0.06 s for GRB 
090510, if r = 10 6 cm. Then A^liv - 0.14 s, and M lC 2 ~ 
9.73 x 10 19 GeV, which is about 8 times of the Planck 
energy -Epianck ~ 10 19 GeV [4j. However, if we increase 
r to 10 7 cm, then At int ~ 0.46 s and At LW ~ -0.26 
s. In this case, high energy photons become superluminal, 
which is against the argument of quantum gravity theories. 
It is a reasonable assumption that ro — 10 6 cm for GRB 
090510, because this is a short burst and its radius should 
be small than long bursts. 

With the above criteria, we try to use the line fitting 
method to predict the LIV effects. In principle, if the lin- 
ear dispersion relation holds, the Atijv/(1 + z) v.s. K (z) 
plot should be a zero-intercept line, whose slope is the 
inverse of quantum gravity energy scale, i.e., 1/Mic 2 . By 
choosing ro of each burst properly, the limits of LIV effects 
for both short and long bursts can be unified respecting 
the linear dispersion relation. The four GRB points can 
be fitted well by one line, if we choose roj = 16.7, 0.1, 28.7 
and 55.0 for GRB 080916c, GRB 090510, GRB 090902b 
and GRB 090926, respectively. In this case, LIV effects 
are calculated in Table [2j and the Atuv/(1 + z) v.s. K(z) 
plot is given in FigfT] The values of ro, which indicate the 
active scale of central engines, are reasonable. The inverse 
of the slope gives M\(? ~ 1.0 x 10 20 GeV, which is roughly 
the same result of the GRB 090510. The energy scale of the 
modified photon dispersion relation is one order of mag- 
nitude higher than the conventional Planck scale. This 
may suggest that the quantum gravity scale may be more 
subtle than one naively thinks, because this quantity is 
model dependent. For instance, the effective quantum en- 
ergy scale depends on the density of _D-particles in the 
D-foam model 0, [1| • It is an interesting future work to 
combine both the quantum gravity model and the GRB 
model together to study the LIV effects. Supposing LIV 
effects are strongly suppressed, which is the assumption 
taken by the mainstream GRB models, the observed time 
delay can well predict the value of ro by using the mag- 
netic jet model, since the tuning of ro is sensitive to the 
fitting. 

4. Discussion and Conclusion 

From Eq.([7]), one can infer that the time lag has ap- 
proximate linear relation with the initial radius (oc r 9 86 ) 
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Table 2: The LIV induced time delay A£liv and quantum gravity energy scale Mic 2 derived from four Fermi-detected GRBs. 
At bs is collected from Ref.0-0- AtLiv = Ai b s — Aii nt , where Atint is calculated by Eq.Q, Eq.([S| and Eq.©. The value of <7o,3 
in each GRB is approximately the bulk Lorentz factor of the jet in unit of 10 3 and is taken as 00,3 ~ 1 [3- r o,7 is chosen as 16.7, 
0.1, 28.7 and 55.0 for GRB 080916c, GRB 090510, GRB 090902b and GRB 090926, respectively 
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Figure 1: The plot of A*liv / v - s - K( z ) f° r f° ur Fermi-detected 
GRBs. 



and the redshift (oc (1 + z) 119 ). It has a weak dependence 
on the photon energy Eq, the peak energy E p , and the 
peak luminosity L >p , which are At- lxx t oc Eq 17 , E p 05 , and 
Lyp 4 , respectively. In addition, the time lag is hardly de- 
pendent on Co. Another model proposed to explain the 
delayed GeV photons was given by Meszaros and Rees. In 
this model, the time lag also depends linearly on tq (See 
Eq.(19) in [2(J). These GRB models also predict the spec- 
trum, so they can be verified by the observed spectrum of 
the whole energy band. 

As mentioned above, in Ellis et al.'s proposal, the linear 
fitting function is written as Ai b s /(l + z) = o-iavK(z) + b. 
In our case, b ~ A^ nt /(1 + z) is not a constant, and we 
can estimate it approximately as 



, n n S J.86 r 0.14 p0.03 pO.16 f-i 



,0.19 



(12) 



which roughly agrees with Ellis et al.'s result b ~ 10 -2 
15]. Since b weakly depends on the redshift, it can be re- 
garded as a distance- independent quantity [17J. However, 
b strongly depends on ro, so the long and the short bursts 
will lead to quite different values of b. If the long bursts 
have roughly the same ro, b can be considered as a con- 
stant, and the linear relation between Af b s /(1 + z) and 
K(z) holds. 

The analysis above gives us the hint that, if we want 



to consider source effects of the GRBs, long bursts with 
small redshifts are preferred. Due to the short cosmolog- 
ical distance, quantum gravity effects can be attenuated. 
With the purpose of enhancing quantum gravity effects, 
short GRB bursts with high redshifts should be selected. 
Therefore, the future observations on short GRBs will im- 
prove the test of LIV effects. If we use an ensemble of 
GRBs with both long and short bursts, the fitting func- 
tion method is not convincing. The intrinsic time delay 
plays an important role. Better knowledge of the intrin- 
sic property of the source will help us to improve the test 
of LIV effects. Inversely, better understanding of quan- 
tum gravity can help us to predict the parameters in GRB 
models. 

In this work, we discussed LIV effects by making use of 
the magnetic jet model. GeV photons are emitted later 
than MeV photons, due to their different optical depths. 
This physical ingredient should be included in probe of 
LIV effects. The neglect of the photon emission mecha- 
nism may lead to misleading results. The constraints of 
LIV effects can be unified for both long and short bursts. 
The calculation of the linear energy dependence of disper- 
sion relation gives M\c 2 ~ 1.0 x 10 20 GeV. Although the 
magnetic jet model itself should be tested by further in- 
vestigations, the analysis of the intrinsic time delay is im- 
portant when we study the photons from the astrophysical 
sources to test LIV effects. 
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